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Abstract 

We have performed an overall fit to the electroweak data with the 
generation blind U (1) extension of the Standard Model. As input data 
for fitting we have included the asymmetry parameters, the particle decay 
widths of Z, neutrino scattering, and atomic parity violation. The QCD 
coupling a s has been constrained to the world average obtained from all 
data except the Z width. On the basis of our fit we have constructed a 
viable gauge model that not only explains R\, and R c but also provides a 
much better overall fit to the data than the Standard Model. Despite its 
phenomenological viability, our model is unfortunately not simple from 
the theoretical viewpoint. Atomic parity violation experiments strongly 
disfavor more aesthetically appealing alternatives that can be grand uni- 
fied. 
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The observation at LEP(T| that the decay widths of the Z to bb and cc 
do not agree with the Standard Model expectations has led to a flurry of 
theoretical activity [§, |], §, §, 0, §]. Various possible explanations have been 
considered. Most of these explanations suffer from at least one defect. Either 
they do not present a complete phenomenologically viable model or they present 
an overall fit that ignores some other experimental data. In this paper we present 
a model that, while aesthetically distasteful, is phenomenologically viable and 
has a much better overall fit to data than the Standard Model. As input data we 
use the various asymmetries and partial widths as measured on the Z resonance 
as well as other data that are constraining. 

First, we describe the philosophy of our model. Couplings of the Z to 
leptons are severely constrained by the current data, so we modify only the 
couplings to quarks. We modify the couplings in a generation independent 
fashion and demonstrate the modifications needed to obtain a good fit to the 
data. This fit is quantified in terms of the total x 2 . We next show that a 
model involving the mixing of the Z with a second more massive boson can be 
constructed. The model is not supersymmetric and requires the existence of new 
quarks to ensure anomaly cancellation. Finally we comment on the constraints 
that the non-observation of such particles and the new gauge boson itself place 
on the model and conclude. 

It is important that any model that purports to explain the problems in 
the bb and cc decay widths of the Z does not introduce problems with other 
processes. Quantities that are measured precisely at the Z are J3J, the mass of the 
Z, the forward-backward asymmetry for leptons (A e FB ), for charm {A C FB ) and for 
bottom (A b FB ) quarks; the asymmetries measured in tau decay (A T and A e ), the 
total width of the Z (T z ), the hadronic production cross section (cx°), the ratio of 
the hadronic to leptonic width (Re), the fraction of the hadronic width that goes 
into charm quarks (R c ) and bottom quarks (Rb)', as well as the left-right beam 
polarization asymmetry (Alr) and left-right forward-backward asymmetries for 
charm (A C (LR)) and bottom quarks (Af,(LR)) |J. In addition there are other 
important pieces of data. The first of these is a s that we constrain to be equal 
to the world average [0] obtained from all data except the Z width; we include 
the measurements from jet counting at the Z |T| since these measurements 
are independent of the couplings of the quarks to the Z itself. Very important 
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are data from lower energy experiments, particularly the measurement of parity 



violation in cesium (Qw) HI! an d thallium (Qw) HH atoms which severely 



constrain the vector couplings of the Z to up and down quarks. The W mass 
{My/) [|rjj severely constrains any shifts in the gauge boson mass spectra and 
finally measurements from neutral current interaction of neutrinos [I5| constrain 
the couplings of up and down quarks to the Z at lower energies. We shall include 
all of these data in our fit. Models that can be favored by the Z data alone are 
disfavored when the rest of the data are included. 

The measured values of R b = 0.2219 ± 0.0017 and R c = 0.1540 ± 0.0074 
deviate by 3.67a and 2.46a, respectively, from the Standard Model predictions. 
The value of R c is 10% lower than the Standard Model value. It is difficult to 
explain these discrepancies by models based on radiative corrections ||. Since 
R q oc g q l + g Q R, where is the left-handed (right-handed) coupling of the 

quark to the Z boson, we need shifts in these couplings due to new physics to 
resolve the Rb and R c anomalies. If the new physics affects only the b and c 
quark couplings, such shifts are difficult to reconcile with the otherwise good 
agreement with the Standard Model for the following reason. Since the QCD 
corrections to the partial decay widths cancel to good accuracy in Rb and R c , a 
shift in Rb and R c changes the total hadronic decay rate into 

r r0 v /, a s {M z ) 2 \ / SR b + 5R c \ 

Thad = n ad x (l + —— + 0(a 2 s ) J * (l + 1 _ r sm_ r sm )> C 1 ) 

where 5R b = R e b xp - R^ M , 5R C = Rf p - R S C M , and T° had denotes the Standard 
Model value of T had before the QCD correction. With a s (M z ) = 0.12, this 
change would shift Tz by —11a from the measured value and, in terms of Ri, 
by —14a. If instead a s (M z ) is extracted by fitting T had to its measured value, 
a s (Mz) would have to be 0.186 ± 0.042, in disagreement with the world average 
of 0.118 ± 0.003. A natural resolution of this problem involves postulating the 
new physics for other quarks too. In particular, if the new physics is generation 
blind, the model is free from the fine tuning problem of flavor changing neutral 
currents, which is a common difficulty in the class of models [Bl |6| that introduce 
new physics only in the heavy generations. 

The simplest way to accommodate these features is to add another U(l) 
factor to SU{2) x U(1) Y of the Standard Model §, 0, g. Mixing between the 
gauge boson X of this extra U(l) with the Z boson of the Standard Model can 
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produce the shifts in the Zqq couplings that are necessary to explain R b , R c and 
a s . The most general generation-blind U(l)x current that is consistent with 
SU(2) x U(1) Y can be written as 

Jx = 9x(q Q Q^Q + quU-fU + q D D^D + qJsfL + q E E-fE + •••), (2) 

where Q and L represent the left-handed quark and lepton doublets, and U, 
D and E are the right-handed up-type quarks, down-type quarks and charged 
leptons, respectively. Summation over generations is understood, and the con- 
tributions from particles other than those of the Standard Model have been 
suppressed. Since U(l)x charges always enter multiplied by gx, we normalize 
them to qQ = — 1 so that five parameters, gx and four charge ratios, specify 
the U(l)x current. The Z — X mixing occurs by Higgs doublets that carry 
U(l)x charges. We assume that there is no higher dimensional Higgs multiplet 
of SU(2). Then the mass matrix may be written as 

M , = ( 2 ) pi 

^ ngzgxv 1 g x V x j 

where g z = #2/ cos6 w with g 2 being the SU(2) coupling and v 2 = (-^IGf)^ 1 if 
the tree level expression for the W mass is unchanged. The gauge eigenstates 
(Z,X) are related to the mass eigenstates (Z M) X M ) by 

Z = Zm cos a — Xm sin a 

( 4 ) 

X = Xm cos a + Zm sin a 

where tan 2a = —2Kg z g x v 2 /(g x V x — g 2 z v 2 /A). The coupling of the Standard 
Model quark of flavor % to the Z gauge boson is given by = gzqi(q^q) with 
q i = (T 3L — Qsin 2 ^^/)^ The lighter mass eigenstate Z M is identified with the 
experimentally observed Z boson. The mixing between Z and X shifts the 
coupling to the observed Z from the Standard Model value by 
$<li — (dx I 'fi , z)o , xi sin a + gj(cosa: — 1), where qxi is the U(l)x charge of quark 
i. If M x ^> Mz, the mixing angle is given by sin a = — ng z v 2 / (g x V x ) and 
therefore 5qi = —Kqxi(v/Vx) 2 - In this approximation there are two parameters 
k, Vx for mixing and four U(l)x charge ratios that are fitted to the data. When 
Mx is comparable with M z , exact diagonalization must be done and the gauge 
coupling gx is included as an independent parameter. 
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Since the Z gauge boson is not a mass eigenstate, the tree-level relation 
M§ = g\v 2 /4 is no longer valid. However, the mass relation Mjy = g 2 v 2 /A is 
not affected. The shift in Mz can be expressed as a shift in the p parameter 
I73fl . Since the Z mass is measured more accurately than the W mass, we use 



the W mass relative to the Z mass as experimental information in comparing 
theoretical predictions with the data. The decrease in Mz is translated into an 
increase in M\y and a decrease in sin 2 9y/. In the large Mx approximation, 



5M W 5M 2 Z COS 9yy 

M w ~ M\ cos 2 6 W - sin 2 6 W ' 



% M 2 Z cos 2 9w- sin 2 ' {:>) 

5 sin 2 9w SM Z cos 2 9w sin 2 9w 

sin 2 9 W ~ M 2 Z cos 2 9 W - sin 2 9 W ' ^ ' 

with 5M 2 Z /M 2 Z = -(M x sina/M z ) 2 . 

Atomic parity violation experiments constrain the vector couplings of the 
up and the down type quarks. For a heavy atom with atomic number Z and 
neutron number N, these experiments measure the charge 
Q w = -2((2Z + N)C lu + (2N + Z)C ld ), where C lq is defined in pf. The 



measured [|T^, |1^] and predicted (TJ, Qvv charges for cesium (Z=55, N=78) 
and thallium (Z=81, N=124) are: 

Q w (Cs) = -71.04 ± 1.81, Q w (Cs) SM = -73.14 (1.16a) 

Q w (Tl) = -114.2 ± 3.8, Q w {Tl) SM = -116.3 (0.55c-). ^ ^ 

Both experiments agree on the sign of the difference between the measured value 
and the Standard Model prediction. These measurements strongly constrain 
any new physics that would further decrease the Qw charge and hence limit the 
values of the U(l)x charges the quarks can have. 

We perform a minimum \ 2 analysis fitting both the shifts in the vector and 
axial couplings of Z and the shifts in My/ and sin 2 0w to the 18 observables 
discussed above. Although the SLD measurement of Alr is inconsistent 
with the LEP measurement, we find no reason to exclude either measurement 
from the fit. Electroweak radiative corrections [15], are incorporated in the 



Standard Model values of these observables. For the Z — X mixing contribu- 
tion to these observables, no radiative corrections are included. Although loop 
corrections can generate kinetic energy mixing between Z and X, such mixing 
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is equivalent to the mass mixing at any fixed q 2 . It makes a small difference 
only when extrapolation is made to different q 2 which is relevant when fitting 
to the low-energy experiments. However, we checked that the Z — X kinetic 
mixing parameters vary by a negligible amount over this range, and so we may 
ignore the extrapolation in the U(l)x current. In the case of the low energy 
parameters, we make only the Standard Model radiative corrections by running 
from q 2 = Mf to q 2 = 0. In computing the electroweak radiative corrections, 
we use a s = 0.118, m t = 175 GeV and l/a(M z ) = 128.75. The new physics re- 
quires a nonminimal Higgs sector. The radiative corrections due to Higgs loops 
are numerically very small. Therefore we approximate the Higgs correction with 
that of the Standard Model by choosing two values (100 GeV and 400 GeV) for 
the Higgs mass. 

In performing the fit we restrict to V x > 750 GeV (see later) and allow the 
leptons to have arbitrary U(l) x charges. We diagonalize the Z — X mass matrix 
exactly. The minimum x 2 is 16 for M H = 100 GeV and the preferred value of the 
U(l)x charges of the leptons is zero. Setting these charges to zero, we have four 
parameters. With fourteen degrees of freedom, our best x 2 is 16 for Mh = 100 
GeV and 15 for Mu = 400 GeV. For comparison, the x 2 f° r the Standard 
Model is 30 for M H = 100 GeV. The U(l) x charges are q v = 2.47 ± 0.33 and 
q D = 1.17±0.46 for M H = 100 GeV, and q v = 2.19±0.69 and q D = 0.905±0.42 
for Mh = 400 GeV. The errors correspond to x 2 = Xmin + 1- See Table 1 for 
the experimental and fitted values of the observables. 

We now build a model based on our analysis. Since the leptons carry no 
U(l)x charge, there are three logical possibilities in constructing a two-doublet 
Higgs model: (1) q v = 2q L - q D ; (2) q v = q Q ; and (3) q D = q Q . We note that 
the fitted U(l)x charges are inconsistent with these possibilities. If the atomic 
parity violation data are excluded, only case (3) is favored by the remaining 
data; the x 2 is 19 (17) for m# = 100 (400) GeV and for 13 degrees of freedom 
(for comparison, the x 2 f° r the Standard Model without these data is 28 for 
16 degrees of freedom). When the atomic parity violation data are included 
the x 2 is 28 for 15 degrees of freedom. As we see no basis for ignoring these 
data, qQ = qo is not feasible and we must introduce three Higgs doublets. 
Anomaly cancellation is challenging and requires many more fermions than in 
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the Standard Model. Q 

A model looks more natural if the U(l)x charge ratios are rational numbers. 
Though this is by no means a requirement, we restrict to this possibility. We 
find that the fitted charges can accommodate such a choice: qQ = — 1, qu = 2 
and qr> = 1. Three Higgs doublets H u , H d and Hi are introduced to give masses 
to the up quarks, down quarks and leptons. Their U{l)x charges are qn u = —3, 
qn d = —2 and qn t = 0. Since the U(l)x charges of the Standard Model fermions 
are not vector-like, the U(l)x gauge symmetry is anomalous and new quarks 
must be added to cancel the anomalies. We add three generations of Standard 
Model-like quarks with opposite U(l)y and U(l)x charges: Q' L = (2, -1/6, 1), 
U R = (1, -2/3, -2), D R = (1, 1/3, -1) under SU(2) x U{1) Y x U(1) X - These 
new quarks, in turn, generate anomalies under SU(2) x U(l)y and their chiral 
partners must be added to make SU{2) x U(l)y vector-like: Q' R = (2, -1/6, 0), 
Ul= (1, -2/3, 0), D L = (1, 1/3, 0). Since all of the added quarks are vector-like 
under SU{2) x U(l)y, there is no contribution to the S parameter |2"Ifl . Since 
twelve quark flavors have been added to cancel anomalies, the QCD coupling is 
no longer asymptotically free. Using the one-loop j3 function, we have checked 
that the coupling remains perturbative up to the Planck scale. 

The new quarks should be heavier than about 200 GeV to avoid detection 
at Fermilab. They can acquire mass through the Higgs doublets: 

\^ l UrE c u + WrDlH? + XsWlDrH^ + KQrUlHu (8) 

where the superscript c denotes charge conjugation. The masses generated by 
these couplings should be of the order of the Standard Model quark masses. 
To make the new quarks heavier, we must introduce additional singlet Higgs 
couplings. These singlet Higgs fields break U(A)x at a scale larger than the 
electroweak scale. Two Higgs singlets (f) and <$' are introduced with the U(l)x 
charges q^ = — 1 and = — 2 so that the new quarks acquire mass through 

+ KUrUl^' + \ 7 ~D R D L <P (9) 

Recently Babu et al || proposed supersymmetric U(l) extension models for Rf, and R c . 
To accommodate supersymmetry qq = qjj was imposed on all models. The good \ 2 that 
they obtained with qQ = qo is mainly due to the neglect of atomic parity violation. If atomic 
parity violation is taken into account, the \ 2 °f their models would be much larger. One of 
their models is specially attractive since it can be embedded into the E(6) grand unified model. 
With atomic parity violation, however, the \ 2 is only slightly better than the Standard Model. 



6 





TVTpa mirpH valnp 


Fit 


Model 


Fit 


Model 






Mh = 100 GeV 


M H = 400 GeV 




2 4963 ± 0032 

uiTi/UU 1 VJ • VJ VJ 


2 500 


2 501 


2 499 

Zj .'-tUU 


2 500 

Zj . W W 


Re 


20 788 ± 032 


20 76 

Zj VJ . ( VJ 


20 78 

Zj VJ . i VJ 


20 76 

Zj . ( VJ 


20 78 

Zj VJ . ( VJ 


u h\ llu ) 


41 488 ± 078 

t: X . t:VJ VJ _1_ VJ • VJ 1 VJ 


41 46 


41 45 


41 45 


41 44 




2210 ± 0017 

VJ . — _ J_ ../ _l VJ . VJ VJ X ( 


2200 


2195 

VJ . Zj X .JtJ 


221 

U.ZZ1U 


2205 

VJ . Zj Zj \J\J 




1540 ± 0074 

UiIUtU —1 — VJ . VJ VJ 1 i 


0.1642 


1649 

VJ . X V Jt: <J 


1626 

VJ . X VJZi VJ 


1 634 

VJ . X V.' <Jt: 




0007 ± 0031 

U.U Jl/ 1 _l_ VJ . VJ VJO X 


1 043 

VJ . Iutu 


1 044 

VJ . X Wtit: 


1 01 7 


1013 

VJ . X VJ X *J 


■^b \- L - , - LL J 


841 ± 053 

VJ . VJ j: X _l — UiUL/U 


9284 

VJ • -J Zj Vjt: 


9297 


9285 


9281 

VJ . JZ VJ X 




0729 ± 0058 

VJ • VJ 1 j-j \J — 1 — VJ . VJ VJ VJ 


0784 

u.u 1 Ut: 


0775 

VJ . VJ \ i '. ' 


0766 

VJ . VJ 1 VJ VJ 


0757 


A (LB) 


606 ± 090 


0.698 


0.690 


0.699 


0.693 


A 


1418 ± 0075 

VJ . X t: X (J _l — VJ . VJ VJ \ 


0.1497 


0.1497 


0.1461 


0.1456 


1 x e 


1390 ± 0089 

V 7 • -L ' 7 V 7 — 1 — V 7 • V 7 V 7 v. 7 ' 7 


0.1497 


0.1497 


0.1461 


0.1456 




0.1551 ±0.0040 


0.1497 


0.1497 


0.1461 


0.1456 


A e 

FB 


0.0172 ±0.0012 


0.0169 


0.0169 


0.0161 


0.0159 


Qw(Cs) 


-71.04 ± 1.81 


-70.74 


-71.78 


-71.03 


-71.36 


Qw(Tl) 


-114.2 ±3.8 


-113.0 


-114.7 


-113.5 


-114.0 


9l 


0.2980 ±0.0044 


0.300 


0.300 


0.3007 


0.3004 




0.0307 ±0.0047 


0.0279 


0.0274 


0.0269 


0.0274 


M w (GeV) 


80.33 ±0.15 


80.43 


80.43 


80.38 


80.37 


x 2 




16 


17 


15 


16 



Table 1: Experimental Jj], |[| [12], [13|, [14], [ijj and fitted values of observables. 
Correlations between the data were included in the fit. Column labeled "Fit" 
shows fitted values of observables for arbitrary U(l)x quark charges (14 d.o.f.). 
Column labeled "Model" gives fitted values for the model discussed in the text 
(16 d.o.f.). All fitted values are for gx = 0.15. The \ 2 f° r the Standard Model 
is 30 (18 d.o.f.) for M H = 100 GeV. 
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When contribution of Eq.([5|) is much larger than that of Eq.(|S|), the new quarks 
are nearly degenerate and a shift in the T parameter JHJ is negligible. 



Since five Higgs fields (H u , H^, Hi, (j) and (j) 1 ) develop vevs, we must ensure 
that they do not result in an unabsorbed Goldstone boson or an axion. We 
introduce self-interactions among the Higgs multiplets to eliminate accidental 
global symmetries that may break down spontaneously. Since two neutral gauge 
bosons of SU (2) x U(1)y x U(1)x absorb two Goldstone modes, we add appropri- 
ate interactions among Higgs fields to give mass to the three remaining modes. 
The following couplings suffice: 

A 8 2 0' c + + \i H d HW c (10) 

The new quarks can be 3 or 3 of SU(3). If (Q', U, D) are assigned to color 
triplets, there is an accidental discrete symmetry, Q — > Q, Q' —> —Q' etc., 
that prevents the lightest new quark from decaying. Then the lightest baryonic 
bound state of the new quarks might be abundant enough to have been detected 
in exotic matter searches p2[ . When they are assigned to color antitriplets, we 



can introduce another scalar singlet <fi and allow the new quarks to decay into Q 
and <p through the coupling QQ'cf). However, the following mass terms are then 
allowed by the gauge symmetries: 

M X Q T L CQ' L + M 2 U T R CU R + M 3 D T R CD R , (11) 

where C is the charge conjugation matrix. These terms result in mixing between 
the Standard Model quarks and the new quarks. They may be forbidden by 
imposing the discrete symmetry mentioned above. We assign an odd parity to 
under this symmetry to maintain the QQ <fi coupling. Since is a singlet 
carrying no U(l)y or U{l)x charge and is stable, it can escape detection in 
terrestrial experiments. The (j) particle could have been produced in the early 
Universe and could contribute to the mass density. The mass and coupling of <fi 
can be adjusted so that it does not overdose the Universe [f23|j . 

We now examine the property of the X boson in our model and some of its 
phenomenological implications. The parameter k in the Z — X mass matrix is 
given by 
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Table 2: 5C/(3) x 5C/(2) x U(1) Y x C/(l)x quantum numbers for matter fields 
in our model. 
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w here (H u ) = v u / y/2, (H d ) = v d / V2 and (Hi) = v t / V2 with 
\l v l + v d + v i = f = 247 GeV. Introducing (0) = V/y/2, (0') = V/y/2 and 
5V 2 = V 2 + 4V' 2 , the parameter V x is given by 5V 2 + 9v 2 + 4v 2 d . Assuming 
vf <C v 2 , the parameters to be fitted are gx, « and V". If the Yukawa couplings 
appearing in Eq.(||) are 0(1), the vevs V and V should be greater than v so 
that the new quarks are heavier than the Standard Model quarks. This implies 
a lower limit V > 250 GeV (or equivalently Vx ^ 750 GeV) that is imposed 
in performing the fit. We restrict gx > 0.1 to guarantee Mx > Mz- Since the 
charges of the up quarks under U(l)x are large, we restrict gx < 0.5 to ensure 
that the coupling strength of the X boson to up quarks (g^+gj^gx/Air = 5^/4^ 
remains perturbative. 

We now discuss a fit to the 18 observables using Mh = 400 GeV. Diago- 
nalizing the mass matrix exactly (which gives the mass of the X boson) in this 
model gives \ 2 — 16, V — 340 GeV, and k = —0.05 for g x = 0.15 fixed. The \ 2 
is not very sensitive to gx- The x 2 increases by 1.5 if gx is varied between 0.1 
and 0.5. With g x = 0.15 fixed, the 90% (95%) C.L. range for V is 250 GeV to 
980 GeV (250 GeV to 1200 GeV) and the 90% (95%) C.L. range for « is -0.01 
to -0.18 (-0.01 to -0.23). The allowed range for k is -1.5 to 1. Therefore k must 
be fine tuned to within 6.8 (8.8) %. 

The mass of the X boson depends on k, V and (almost linearly) on gx- For 
small k, M 2 X ~ g 2 x (5V 2 + 6v 2 ). The 90% (95%) C.L. range for M x is from 125 
GeV to 341 GeV (125 GeV to 412 GeV) for g x =0.15 and M H = 400 GeV. For 
g x = 0.5 and M H = 100 GeV, the 90% (95%) C.L. range for M x is from 420 
GeV to 820 GeV (420 GeV to 1070 GeV). The X boson can be produced in pp 
collisions at the Tevatron and detected in the dijet final state. For gx = 0.15, 
the expected production rate is considerably below the limit set by the CDF 
group |24| for all values of V. For g x = 0.5, the values of M x & 750 GeV 
(V& 610 GeV) are excluded. For g x = 0.3, the region 320 GeV <>M X ^ 520 
GeV is excluded. 

To summarize, on the basis of an overall fit to all electroweak data we have 
built a viable U(l) extension of the Standard Model. While the fit to data has 
been greatly improved, the model lacks aesthetic appeal. The X boson may be 
accessible by the experiments at Fermilab in the future. 
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